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Molecular electronics on silicon has distinct advantages over its metallic counterpart. We describe 
a theoretical formalism for transport through semiconductor-molecule heterostructures, combining 
a semi-empirical treatment of the bulk silicon bandstructure with a first-principles description of 
the molecular chemistry and its bonding with silicon. Using this method, we demonstrate that the 
presence of a semiconducting band-edge can lead to a novel molecular resonant tunneling diode 
(RTD) that shows negative differential resistance (NDR) when the molecular levels are driven by 
an STM potential into the semiconducting band-gap. The peaks appear for positive bias on a p- 
doped and negative for an n-doped substrate. Charging in these devices is compromised by the 
RTD action, allowing possible identification of several molecular highest occupied (HOMO) and 
lowest unoccupied (LUMO) levels. Recent experiments by Hersam et al. pj support our theoretical 
predictions. 
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Traditionally molecular electronic efforts, both theo- 
retical 0, and experimental 0, EJ 0, 0, HI > have been 
driven by thiol-gold chemistry to molecules bonded to 
gold substrates. However several recent experiments 
have demonstrated the feasibility of attaching various 
molecules on clean silicon surfaces The development 
of molecular electronics on silicon is particularly impor- 
tant for two reasons. Firstly, it will enable the develop- 
ment of integrated devices that can utilize the powerful 
infrastructure provided by the silicon-based IC industry. 
Secondly, unlike gold, silicon has a bandgap that one can 
take advantage of to design a new class of resonant tun- 
neling devices with pos sible applications in logic ^1 an d 
low-power memory 111. In view of these significant po- 
tential payoffs we believe it is worthwhile at this time to 
develop models that can be used to analyze the electrical 
characteristics of such silicon-based molecular devices. 

The purpose of this paper is to present (i) a general 
formulation suitable for modeling silicon-based molecu- 
lar devices; (ii) realistic principles for designing RTDs 
based on such structures; and (iii) showing ways to map 
out the molecular energy spectrum that is not realizable 
in normal experiments with gold contacts. Molecular de- 
vices on silicon require a formulation that can account for 
the bandgap, surface band-bending and surface recon- 
struction. One problem commonly encountered is that 
standard quantum chemical basis sets describe molecu- 
lar energy levels well but provide a very poor descrip- 
tion of the semiconductor bandstructure. We present a 
scheme that can be used to integrate the two distinct 
systems seamlessly and present results that combine ab 
initio treatments of the molecule with a semi-empirical 
description of the silicon bands. The same approach can 
be used to integrate the molecule with more advanced 
treatments of the silicon substrate. 

Our results show that we can expect to see NDR using 
molecules like styrene or TEMPO (2,2',6,6'-tetramethyl- 
1-piperidynyloxy) on p-silicon at realistic positive sub- 
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FIG. 1: Schematic description of molecular RTD involving 
a silicon band-edge and an STM tip. The extended device 
Hamiltonian H includes the molecule and a part of the sili- 
con surface relevant to the bonding with the molecule. The 
influence of bulk silicon and STM contacts are incorporated 
through self-energy matrices Ei,2. Bottom: (a) Equilib- 
rium band-alignment in p-silicon-molecule-metal heterostruc- 
ture corresponds to a flat Fermi energy Ef (dashed line) near 
the silicon valence band-edge (VB). Shown also are the sili- 
con conduction band (CB), the metal band (right), and the 
molecular energy levels, (b) For positive substrate bias, the 
levels move up until the HOMO levels leave the bulk silicon 
valence band into the band-gap, leading to a sudden drop in 
conductance and a corresponding NDR in the I-V. 



strate voltages, when current flow occurs through the 
HOMO levels. A schematic description is shown in Fig.^ 
Recent experimental observations are in good agree- 
ment with these results. By contrast, it appears that 
much larger negative voltages will be needed to observe 
these effects on n-type silicon where conduction takes 
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place via the LUMO levels. This is an important observa- 
tion that deserves careful attention, since little work has 
been reported to date for electronic conduction through 
the LUMO levels. 

Formalism. We will use a method by which a bulk sil- 
icon bandstructure calculation (e.g. effective mass 01 
or tight-binding sp 3 s* [13|) can be connected up with a 
detailed ab-initio treatment of the molecule, without in- 
troducing spurious basis-transformation related effects at 
their interface. We accomplish this connection through 
the real-space Green's function obtained on a surface sili- 
con atom, which can be looked upon both as a scooped- out 
part of the silicon crystal, as well as an extended part of 
the molecular device. The connection is established as: 



g(E) = [ESo-Ho-^E)}- 



(1) 



where So and Hq are the ab-initio overlap and Hamil- 
tonian matrices for a surface silicon atom, Si is the de- 
sired self-energy coupling it to the bulk, and g is the 
silicon surface Green's function. We start with an ap- 
propriate semi-empirical description of the bulk silicon 
bandstructure, and use a real-space, recursive formalism 
along the (100) direction to get the semi-empirical surface 
Green's function 0. We then use Eq. \i\to adjust the 
self-energy Si to achieve the best match between the sur- 
face Green's functions (ab-initio and semi-empirical) in 
real-space (Fig.|2J, including if required a suitable recon- 
structed cluster for unpassivated surfaces. We thus have 
an ab-initio description of the molecule and its bonding 
with the silicon surface that is consistent with a semi- 
empirical description of the bulk and surface properties 
of the silicon crystal . 

We use a recently developed method 0,0] to calculate 
the current, coupling the nonequilibrium Green's func- 
tion (NEGF) formulation of transport ^| with an LDA 
molecular Hamiltonian calculated using GAUSS1AN98 
in a 6-31g* basis: 



G = [ES - H - Si - S 2 ] 
T = Trace [r 1 GT 2 G + ] , r\ 2 



Si.: 



J l,2 



{2e/h) / dET(E)[f 1 (E)~f 2 (E)} 
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where S, H are the molecular overlap and Hamilto- 
nian matrices and fx,2 are the contact Fermi ener- 
gies. The calculated current includes contributions 
both from thermionic emission over the vacuum bar- 
rier, as well as quantum mechanical tunneling through 
it. We include the effects of an incoherent back- 
ground inside the molecule through an additional phase- 
randomizing voltage-probe |2fj| , incorporated in the gen- 
eralized Landauer-Biittiker formula |2Jj . 

It is well known |22j that unpassivated silicon has sur- 
face states in the band-gap, as we indeed observe (Fig.[2J). 
Some of these states are eliminated once the molecule 



bonds to the surface silicon dimers |23j. The surface 
states in unpassivated silicon would influence the STM 
current only if they can be replenished from the bulk 
through inelastic processes occuring in the contact. We 
ignore such effects in this paper. 
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FIG. 2: Results for Si(100) density of states (DOS) in a hybrid 
basis: (a) bulk DOS in (I) effective mass theory and (II) ex- 
tended Huckel theory (EHT) with parameters fitted to match 
the bulk silicon bandstructure [T3 |. The surface DOS match 
in energy domain for (III) EHT and (IV) 6-31g* bases, (b) 
The surface DOS match in real space as well, plotted here at 
the Fermi energy, ao denoting the Bohr radius. 

In any STM measurement on a molecule, it is impor- 
tant to deconvolve the tip and the molecular density of 
states (DOS) so that one can unambiguously identify fea- 
tures intrinsic to the molecule. Modeling tip-sample in- 
teractions requires a precise knowledge of the experimen- 
tal tip structure, as well as its exact position over the 
molecule. It is generally believed [24| that the tunneling 
current drawn by a tungsten tip is dominated by the d- 
electrons of loose tungsten atoms at its end. We model 
the tip-sample interaction by evaluating the largest cou- 
pling matrix element between the styrene HOMO wave- 
function and the d-electrons of a single tungsten atom 
sitting at varying distances from it. The coupling is ob- 
tained using LDA in a LANL2DZ basis, and yields a 
broadening ~ 300 meV when the STM almost touches the 
molecule and 12 meV when the tip sits about 2 A away. 
The STM self-energy S 2 is obtained using Fermi's Golden 
rule, with a tip DOS - 0.25/eV. We find that the cou- 
pling can vary considerably depending on the lateral po- 
sition of the tip atom relative to the molecule, due to ac- 
cidental symmetries in the overlap matrix elements. This 
needs further consideration both from theory and exper- 
iment, given that the current level is influenced substan- 
tially by the STM tip. 

Much of molecular RTD action hinges on the bias- 
dependence of the molecular levels and broadenings. Us- 
ing the semiconductor bulk band-edges as reference, the 
electronic energy levels and wavefunctions, the silicon 
band-bending, as well as the tip electrochemical potential 
vary under bias (Fig. ^) . Under bias, the molecular levels 
move relative to the silicon band-edge, because (1) a large 
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part of the applied bias drops across the tip-to-substrate 
gap due to the large dielectric constant of silicon relative 
to the molecule and the vacuum gap; and (2) the field 
inside the molecule is poorly screened by the absence of 
metal-induced gap states (MIGS) that would normally 
exist with gold contacts (this creates a low quantum ca- 
pacitance which causes a sizeable potential drop within 
the molecule HSIH^)- At a critical voltage determined 
by the electrostatics, the HOMO level in Fig. ^ reaches 
the bulk silicon valence band-edge (VB) after which it en- 
ters the Si band-gap and its transmission drops abruptly, 
leading to a prominent NDR in the I-V of the molecule. 
Thus the self-consistent potential U sc { is the most impor- 
tant part of the modeling. 

The equilibrium band-bending in silicon is given es- 
sentially by the workfunction difference between silicon 
and tungsten, and is distributed over a depletion width 
~ 1 nm. We obtain the volta ge-d ependence of the silicon 
band-bending using MEDICI 27], which agrees with sim- 
ple estimates based on the depletion-width approxima- 
tion 1281. a s well as surface photovoltage (SPV) measure- 
ments 29]. The band-bending affects the silicon surface 
Green's function through the onsite energies of the silicon 
atoms distributed over a depletion width. Band-bending 
turns out however, to be relatively unimportant, since it 
is transparent to electron tunneling, so that the NDR is 
determined by the bulk silicon band-edges. The Laplace 
part of the potential drop is obtained using MEDICI [22j 
by simulating a silicon surface separated from tungsten 
by a vacuum gap and a molecular dielectric with a dielec- 
tric constant ~ 2. Stark effects in the molecular DOS are 
included in GAUSSIAN98 by applying a molecular field 
extracted from the Laplace potential. 

Results. Eq. ^can be used to obtain the coupling of a 
molecule to silicon given a bonding geometry and a sur- 
face structure. The actual bonding geometries vary [30| . 
and are sometimes controversial pU| . so we will adopt 
specific models in this paper for illustrative purposes. We 
start by simulating transport through styrene grown on a 
p-doped (N A = 5 x 10 19 cm -3 ) H-passivated Si(100) sur- 
face (bonding shown in Fig. 0Jd). Our LDA energy levels 
and wavefunctions of styrene agree well with previously 
published first-principles results for the above bonding 
geometry |32j. We use our hybrid basis scheme to con- 
nect this geometry to bulk silicon, described here within 
an effective mass model. The effective mass model is 
used for simplicity and describes the band-edge proper- 
ties correctly (Fig El) where most of our physics lies; fu- 
ture publications |lfil ] will generalize this to more sophis- 
ticated models for bulk silicon. Fig. |21 shows a colorplot 
of the transmission of styrene on p-Si under bias. For 
increasing positive substrate voltage the molecular levels 
increase in energy until a HOMO level reaches the silicon 
valence band-edge and its transmission gets cut-off there- 
after. This leads to an NDR in the I-V (Fig.g^), calcu- 
lated for an STM tip almost touching the molecule. The 
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FIG. 3: Calculated transmission (colorplot) through styrene 
absorbed on p-silicon. The silicon is hydrogen-passivated and 
doped at Na — 5x 10 19 cm -3 , giving a Fermi energy Ef ~ 0.1 
eV below its valence bandedge (Ev)- Styrene HOMO peaks 
shift under bias, and their transmission disappears when they 
reach the band-edge Ey. The transmission is plotted for three 
candidate voltage values, marked by crosses in the colorplot. 

reduced transmission also leads to an overall current level 
that is much smaller for positive substrate bias. Multi- 
ple HOMO levels contribute to multiple NDR peaks. For 
negative bias, there is no such molecular NDR, since the 
HOMO levels move deeper into the valence band, and the 
LUMO levels lie outside the bias window. An increased 
tip-to-sample vacuum gap (2 A) causes a smaller voltage 
drop across the molecule, postponing the onset of the 
NDR peaks (dashed line). The NDR peak position from 
a level Eq is roughly given by {Ey — E$)/ri, where r\ de- 
notes the average molecular potential under unit applied 
bias (ij ~ 0.5 for an STM touching the molecule, decreas- 
ing approximately linearly with increasing vacuum gap) . 

For molecules on n-silicon, NDR is expected at nega- 
tive substrate voltages, mediated through the LUMO lev- 
els. The LDA LUMO levels for styrene and TEMPO are 
a few volts higher than Ec, pushing the first NDR peak 
outside the usual bias range. Fig.^J; shows the results for 
n-styrene, with the LUMO levels lowered artificially by 
about 1.5 eV. Electron addition levels are usually contro- 
versial, and could differ from theoretical predictions due 
to various effects associated with charging, correlation or 
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FIG. 4: (a) Calculated I-V characteristics for styrene on p- 
doped H-passivated Si(100) with an STM almost touching 
the molecule. Increasing the air gap to 2 A (dashed line) 
changes the overall current and postpones the NDR action; 
(b) the bonding geometry [§3 between an end carbon atom 
on styrene and a c(4x2) Si(100) buckled dimer; I-Vs for (c) 
styrene on n-doped H-passivated Si(100) and (d) TEMPO on 
n-doped clean Si(100), N D = 5 x 10 19 cm" 3 . 



image effects. The theoretical NDRs compare well with 
room temperature STM measurements by Hersam et al. 
Q on TEMPO and styrene, both single molecules and 
monolayers. The experiments are performed on clean 
Si(100) surfaces, with a different surface reconstruction 
(2x1), bonding mechanism (cycloaddition) , and band- 
bending |33(, but the corresponding theoretical results 
are qualitatively similar (Fig.^Jl). 

The origin of NDR in our simulations is qualitatively 
different from the NDR mediated by tunneling through 
localized atomic states arising from a sharp structure in 
the tip DOS j3^|. NDR due to a tip to level resonance 
would show up either for both bias directions for both 
p and n-type samples (if the tip structure interacts with 
both HOMO and LUMO levels), or for one bias direc- 
tion for only one of the samples (if it interacts with just 
one of them). The consistent observation of a one-sided 
NDR, at positive substrate bias in p-Si and negative in 
n-Si yj/, can only be explained by interactions between the 
molecular levels and the silicon band- edge. 

The multiple NDR peaks arise from the multiplicity of 
levels sequentially crossing the relevant band-edge. This 
is a significant improvement over gold-based experiments, 
where large molecular charging energies (~ 2 eV) prevent 
an STM electrochemical potential from crossing more 
than one level. So far, experimental molecular conduc- 
tances exhibit very few (typically zero Q , one 0, @ or 
two 0) peaks, in contrast to systems with smaller charg- 



ing energies |35|. In molecular RTDs, however, the prox- 
imity of the band-edge to the Fermi energy dramatically 
reduces this charging effect, because any charging due 
to electron/hole addition to a level by the stronger cou- 
pled substrate is counterbalanced almost instantly by its 
removal due to the dramatic reduction in substrate cou- 
pling at the bandedge. The low charging allows the obser- 
vation of multiple NDR peaks that can be correlated with 
the molecular levels. Indeed, we can use a combination 
of p and n-type experiments on a single molecule to map 
out its HOMO-LUMO spectrum. For instance, the ex- 
perimental peak voltages A„ iP and Fermi energies Ef n j p 
for n and p-doped substrates give a HOMO-LUMO gap 
« rj (A„ + A p ) + (1 - rj) (E fn - Ef p ), including equi- 
librium tip-substrate band-alignment effects. Charging 
can also lead to a hysteresis in the I-V, associated with 
the changeover of the substrate coupling from strong to 
weak at the NDR position. Such hysteresis effects are 
observed in RTDs |36j, but will also be compromised by 
the low charging in these molecules. For realistic device 
parameters, the predicted hysteresis is small, and could 
be washed out by scattering. 
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